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Abstract 

In this paper, we investigate the performance of cooperative spectrum sensing (CSS) with multiple 
antenna nodes over composite and generalized fading channels. We model the probability density 
function (PDF) of the signal-to-noise ratio (SNR) using the mixture gamma (MG) distribution. We 
then derive a generalized closed-form expression for the probability of energy detection, which can 
be used efficiently for generalized multipath as well as composite (multipath and shadowing) fading 
channels. The composite effect of fading and shadowing scenarios in CSS is mitigated by applying an 
optimal fusion rule that minimizes the total error rate (TER), where the optimal number of nodes is 
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derived under the Bayesian criterion, assuming erroneous feedback channels. For imperfect feedback 
channels, we demonstrate the existence of a TER floor as the number of antennas of the CR nodes 
increases. Accordingly, we derive the optimal rule for the number of antennas that minimizes the TER. 
Numerical and Monte-Carlo simulations are presented to corroborate the analytical results and to provide 
illustrative performance comparisons between different composite fading channels. 

Index Terms 

Energy detection, multipath/composite generalized fading channels, cooperative spectrum sensing, 
square law selection. 


I. Introduction 

The need for efficient utilization of spectrum has become a fundamental requirement in modern 
wireless networks, which is mainly due to spectrum scarcity and the ever-increasing demand 
for higher data rate applications and Internet services [1]. Cognitive radio (CR) networks is 
a particularly interesting proposal that has been developed to mitigate the spectrum scarcity 
issue. CR networks can adapt their transmission parameters according to the environment [2]. 
Cognitive radios have been shown to be very efficient in maximizing spectrum utilization due to 
their inherent spectrum sensing capability. In a CR network environment, users are categorized as 
either primary users (PUs) or secondary users (SUs). PUs are the ones who have been assigned 
spectrum slots, and hence, have higher priority. On the other hand, SUs are opportunistically 
accessing the under-utilized spectrum. Three detection techniques are commonly used for spec¬ 
trum sensing in CR networks, namely energy detection (ED), matched filters, and cyclostationary 
detection [3]. In energy detection, which is the most common detection techniques, due its 
low implementation complexity, the presence of a PU signal is simply detected by comparing 
the output of the energy detector with a certain threshold that depends on the SNR [4]. In 
matched filters, exact information about the transmitted signal waveform, such as its bandwidth 
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and modulation type, is required. Cyclostationary detection, on the other hand, uses statistical 
properties of the transmitted signals to enhance the probability of PU detection. 

Recently, the ED performance over small scale fading channels, such as Rayleigh, Rician, and 
Nakagami-m is studied in [5], [6]. In addition, unified ED performance analysis over generalized 
multipath fading channels has been studied by the utilization of the k — /i and k — // extreme 
fading channels, which outperforms the characterization capabilities of Rician, Nakagami-m 
and Rayleigh distributions [7]. However, in addition to small scale fading, in most scenarios, 
the received signal is also degraded by large scale fading (shadowing effect). Therefore, there 
is a need to analyze the CR performance over composite fading channels, which take into 
account large scale and small scale fading. It is worth noting that a direct analysis of the ED 
performance in composite fading channels is rather tedious, since composite fading models can 
only be represented by infinite integrals. The probability of detection of the ED is derived 
over Nakagami-lognormal (NL) fading channels in [13]. However, the offered solution is not 
represented in closed form and the impact of fading and shadowing effects is investigated 
numerically. Several alternatives that characterize composite fading channels are presented to 
provide simplified performance analysis for the CR networks. For example, in [8]-[12], the /C 
distribution is utilized to study the ED performance over RL channels. 

One the most effective techniques to mitigate the effect of fading is cooperative spectrum 
sensing (CSS), where N users sense the spectrum independently and send their soft or hard 
decisions through an imperfect feedback channel to a fusion center (FC) for a global decision 
[14]. In CSS, there are two main decision combining schemes, namely soft decision combining, 
where the SUs send their local observations to the FC, and hard decision combining, where 
the SUs send their local decisions to the FC. This paper considers the hard decision combining 
scheme. In [15], CSS optimization over the Rayleigh fading channel has been studied, assuming 
a perfect feedback channel. In [16], Quan et. al. studied the performance of CSS under the 
Bayesian criterion, which accounts for the costs of probabilities of miss-detection and false alarm. 
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However, the feedback channel is assumed to be perfect and the study is limited to Rayleigh 
and Suzuki fading channels. In [17], Lee studied the performance of CSS over Rayleigh fading 
assuming imperfect feedback channels and hard decision combining. Recently, CSS over n — fi 
fading channels has been studied in [18]. However, this model does not represent composite 
fading channels, and the authors considered only the OR fusion rule with a perfect feedback 
channel. 

Although CSS has been investigated extensively in the literature, none of the reported results 
has considered CSS optimization over composite and generalized fading channels. In this paper, 
we propose a mixture gamma (MG) based approach to derive a new closed form expression of 
probability of detection over generalized fading channels. The MG model [19], [31] has been 
proposed as an alternative to various non-composite and composite fading channels: Lognormal, 
Weibull, Rayleigh/Lognormal (RL), NL, /C, KLq, 77— /z, k—/ 1 , Hoyt, and Rician channels. The MG 
model is both accurate and flexible to represent all aforementioned fading channels. It is worth 
noting that the probability of detection of the MG model has been derived in [19]. However, the 
solution provided therein is limited to an integer /3 k, which is a parameter in the MG model (4). In 
this paper, we provide a different approach and derive the probability of detection for both integer 
(3k and arbitrary real 3k- In addition, the derived probability of detection is extended to show the 
effect of diversity reception by employing the square-law selection technique (SLS). Moreover, 
we derive an optimal fusion rule where the optimal number of nodes is derived analytically 
under Bayesian criterion. Furthermore, we found a TER floor, whereby increasing the number 
of antennas after this floor will not reduce the TER. As such, we derive the optimal number 
of antennas that achieves the optimal performance. We stress out that results of this paper are 
applicable to all fading channels. Our derivations incorporates both cooperative sensing, diversity 
reception, and assumes erroneous feedback channels. More precisely, our main contributions are 
summarized as follows: 

• The probability of detection of the energy detector using the MG model is derived with and 
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1 M 



Fig. 1: Cooperative spectrum sensing network configuration. 


without diversity reception. 

• An optimal fusion rule is derived under Bayesian criterion for CSS optimization over 
generalized fading channels. 

• A closed form expression of the optimal number of antennas M* is derived. 

The rest of the paper is organized as follows. Section II gives a detailed description of the 
system model, while in Section III, we derive the probability of detection using the MG model. 
In Section IV, we extend our derivation to diversity reception using the SLS technique. Optimal 
fusion rule and optimal number of antennas over generalized fading channels are presented in 
Section V, followed by conclusions in Section VI. 
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II. System Model 


We consider a cooperative spectrum sensing network configuration, shown in Fig. 1. The 
configuration consists of one primary user (PU) and N secondary user (SU) nodes, each equipped 
with M antennas. The SU nodes report back their hard decisions via imperfect feedback channels 
to a fusion center (FC), which may represent a base station or a cell site. We assume that the 
channel gains are independent and identically distributed ( i.i.d ) and are modeled using the 
generalized MG fading model, where % = 1,.., N, and j = 1,..., M. 

The received signal copies at the i th SU node and j th antenna can have two possible hypothe¬ 
ses, modeled as 

( 1 ) 

Ho ■ yfit) = hf{t)s{t ) + vf(t ), 


where H 0 and II \ represent the absence and presence of the PU, respectively. s(t) corresponds to 
the transmitted signal from the PU, with energy E[|s(t)| 2 ]= E s . Vj‘\t ) ~ CM((), a'l) is circularly 
symmetrical complex additive white gaussian noise (AWGN). 

Each i th SU node uses an energy detector for the amplitudes \y 'j\f = \ of the received signal 
and compares it to a threshold A. Therefore, the output of this process for each antenna can be 
written as 


Z f = 



2 Hi 

Ho 


A, 


( 2 ) 


where we have dropped the time index t for simplicity. 

For this set-up, the CSS scheme performs following steps: 

• Each i th SU calculates the decision statistic for all M antennas and uses the SLS technique 
to decide on the presence or absence of a PU. 

• The binary hard decision of each SU is sent to the FC via an erroneous feedback channel. 

• The FC applies a fusion rule to the binary decisions received from all SUs and a final 
decision is made. 
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The conditional probabilities of detection and false alarm are calculated using [4] 

Pd = Qu(V^T>V\i) 


Pf = 


r(«,¥) 


( 3 ) 


/ i» 

where u is time-bandwidth product, Q u (-,-) is the generalized Marcum-Q function [24], T(-, •) 
is the upper incomplete gamma function [25, eq. 8.35], T(-) is the standard gamma function [25, 


eq. 8.31], X n = A/ay 2 is the normalized threshold, and 7 ] 1 = 
of the i th PU-SU link. 


(0 _ \ h f I 


2 ai 


is the instantaneous SNR 


Since the probability of false alarm is based on the null hypothesis, it remains the same 
regardless of the fading scenario. Thus, in the subsequent sections, we focus on the derivation 
of the probability of detection for the MG distribution with and without diversity. 


III. Probability of Detection over Composite Fading Channels with No 

Diversity Reception 

The MG distribution has been shown to accurately represent several composite and non¬ 
composite fading models, and it is written as [19] 


c 


, . ak , x . fa 1 C k%\ 

W = V —(—) exp(-), 

^ 7o 7o 7o 


( 4 ) 


where C denotes the number of mixture components, ak is the mixing coefficient of the k lh 
component having the constraints a k > 0, with a k = 1- 7o is the avarage SNR. /3 k and 
( k correspond to the scale and shape parameters of the k th component, respectively. Using the 


definition, the probability of detection over MG model is given by 

c 

Pd ,mg 


[ Qui^Vxn-f^e-^dx. 


( 5 ) 


Following lemma 1 in [33, eq. 8], and after some simplifications, (5) is solved in exact infinite 
series representation as follows 


„ _^a*.^(v^V‘r(ft + i)r(u + i,i ) 1 

Pl-.MG - 2 _^— \ 


k= 1 i = 0 l\T(u + /)(! + “) 


C k\ 0k+l 


( 6 ) 
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Although, the above solution is expressed as an infinite series, it converges after relatively few 
terms. However, for the sake of completeness, it is essential to determine a corresponding simple 
and accurate truncation error of ( 6 ). Following lemma 2 in [33], a closed form upper bound for 
the truncation error, e t , can be derived as 


c 


^2a k T((3 k ) ^(v/2) 2^(0* + l)T(u + l, £ 

£t ~ 2^ ~ 


k= 


I 7 » (i) & U «r( M + ;)(2 + 2£)‘ + ' 


(7) 


' 70 / ^ \ 70 ' 

Additionally, an exact finite series representation is derived given that /3 k is limited to integer 

values. With the aid of theorem 1 in [33, eq.3], and after some simplifications, (5) is derived as 

c 


-p.int 

r d,MG 


Ofc ,-r(/3fc)r(w, 

^7o 


k =1 


'70 ' 


fc-r (Xrm) 1 F 1 (l + l,u + 1 >T ^ ET ) 
- 

M !(^) /3fc " 1 2“(l + ^) Z+1 ef 

V 70 ' ' A'n / 


4+(4—) 
v 70 y 


( 8 ) 


70 ‘ 


where the function 1 Fi(.,.,.) is the confluent hypergeometric function [23, eq. 9.210.1]. 

In Fig. 2, the complementary receiver operating characteristic (ROC) curves over different non¬ 
composite and composite fading channel is plotted using ( 6 ). The term ’Numerical Integration’ 
refers to numerically integrating (5) via the trapezoidal integration routine in MATLAB. It 
can be observed that there is a perfect match between the analytical form derived in ( 6 ) and 
its corresponding simulation, where as shown, the curves span a wide average SNR range. 
The addressed scenarios were all approximated by two mixture components, 0 = 2. The 
approximated fading channels and their corresponding MG coefficients can be found in Table I 


[30], 


IV. Probability of Detection over Composite Fading channels with Diversity 

Reception 

In the SLS scheme, for each i th SU node, the branch with the maximum 7 j is selected as 
follows [5] 


7 SLS = max(7 j),j = 1,.., M , 

j 


(9) 
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TABLE I: MG coefficients for different composite and non composite fading channels 


Channel 

Param 

a 

P 

VC 

MSE 

Lognormal 

C = 1 

0.8795306 

0.1204694 

21.47391625 

19.29669267 

0.04634621 

0.06526572 

2.18E-04 

Weibull 

m — 4 

0.4163066 

0.5836934 

2.547627 

6.410809 

0.276866 

0.1887993 

8.14E-05 

Rayleigh 

- 

0.2803494 

0.7196506 

0.9124631 

1.3046339 

0.4795083 

0.9339172 

- 

Nakagami-m 

- 

77i m 

r(m) 

m 

J_ 

771 

- 

Ray/logn 

f = 1 

0.2889985 

0.7110015 

0.9667361 

0.972492 

3.908869 

1.62796 

1.18E-06 

Ray/logn 

C = 1-5 

0.7229848 

0.2770152 

0.9223047 

0.8499458 

0.8351385 

3.0324419 

- 

Ray/logn 

c = 0.5 

0.3491298 

0.6508702 

0.9117919 

1.225578 

0.6304104 

1.010205 

1.43E-06 

Nak/logn 

m = 2 

C = 0.5 

0.6569638 

0.3430362 

1.9707105 

2.5034565 

0.4273802 

0.5267421 

1.55E-06 

Nak/logn 

m = 4 

C = 0.5 

0.7775037 

0.2224963 

4.0356456 

5.2938705 

0.2247047 

0.2556295 

3.76E-06 


where M denotes the number of antennas at each SU node. Under H 0 , the probability of false 
alarm for the SLS scheme can be expressed as 

Pf,SLS = 1 — P r (7SXS < ^n\H 0 ). (10) 

Substituting (9) in (10), we obtain 

Pf,SLS = 1 - Pr(max(7i,7 2 , ..,7m) < K\H 0 ). (11) 

Using order statistics [29], P/,sls is written as 

Pf,SLS = 1 - [1 - Pf] M - (12) 
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Fig. 2: Complimentary ROC for single SU over various fading channels. 


Similarly, the unconditional probability of detection over AWGN channel is obtained by 

M 

Pd,SLS = 1 — ]^[ [l — Qu{\J2jj, \/~K) • 


(13) 


3 = 1 

Averaging (13) over (6), we obtain the unconditional probability of detection under the SLS 
scheme, P<i,sls , expressed as 


M 

Pd,SLS = 1 — [1 ~ Pd,M g] ■ 

3 = 1 


(14) 
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V. Optimal fusion rule and optimal number of antennas in Composite Fading 

Channels 

Cooperative spectrum sensing can highly decrease the probability of miss-detection in fading 
channels. In this section, the performance of cooperative spectrum sensing is studied over various 
fading channels. Practically, the channel between the SUs and the FC is imperfect. Therefore, the 
hard decisions sent by each SU node are affected by the quality of the feedback channel. Here, 
a binary symmetric channel with error probability q is assumed. The probabilities of detection 
and false alarm over erroneous feedback channels are given by [17] 


p} = (i — q)Pi + 9(1 - A), (15) 

P'd,MG = (1 ~ + (9)(1 - Pd,MG)- (16) 

For the SLS diversity reception, we substitute (12) and (14) in (15) and (16) respectively, yielding 

P',,SLS = (1 - (1 - A)")( 1 - 9) + (1 - Pf) M 1, (17) 

and 


P d,SLS = 

M M (18) 

(1 — g) (1 — n (1 — Pd, mg)) + q n (1 — Pd,Mc)■ 

3 = 1 i =1 

Assuming i.i.d diversity branches, (18) can be expressed as 


P'd.SLS = (1 - 9) (1 - (1 - Pd,MG) M ) + 9(1 - Pd,Mc) M ■ (19) 

The probability of false alarm Qj- and the probability of detection Q d at the FC using /c-out-of-A r 
rule is thus given by 


N 


Qf - 


n=k 


^ N ^ 


V n ) 


( P'f,SLs ) n (1 — ( P'f,SLS )) 


N-n 


— 1 — B F (k — 1, N, P^sls) i 


( 20 ) 
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and 


Q d — t: 


N { N \ 


( P'd,SLs ) n (1 — ( P'd,SLS )) 


N—n 


n=k y 77/ j 

= 1 — B f (/c — 1, iV, P'd SLs)i 


( 21 ) 


where B F is the binomial cumulative distribution function [26]. Since P m = 1 — Pd, the 
probobility of miss-detection at the FC is given by 


Qm ^ ' 


N ( N \ 


(1 P'd ) SLs) n (P'd ) SLS ) 


N—n 


n=k y Tl j 

= B F (fc — 1, N, 1 — P'd SLS )• 


( 22 ) 


A. Optimal Fusion Rule under Bayesian Criterion 

The performance of CSS in composite and non-composite fading environments can be im¬ 
proved by applying an optimal fusion rule to decrease the TER. The problem is formulated under 
the Bayesian criterion which accounts for the cost of miss-detection and false alarm. Such cost 
can be set by a regulatory or based on past experiences which makes the preceding analysis both 
practical and widely applicable to any fading channel. We consider the problem of minimizing 
the Bayesian risk of CSS, which is given by [16] 

min R(k ) = W f Q f + W m Q m , (23) 

1 <k<N J 1 

where W m and Wf are the cost of miss-detection and false alarm, respectively. From (20) and 
(22), the Bayesian risk R(k ) can be written as 

m = Wf~ WfB F (k - 1, N, P' LSLS ) 

+ WmBp^k — 1, N, P' d SLS ), (24) 

Let R 0 (k) be defined as 

R 0 (k) = W f B F (k - 1, N, P f ) - W m B F (k - 1, N, P d ), (25) 
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then we have 


d R(k) 
dk 


^ N ^ 


\Wm(P'd,SLs) k ( 1 P'd,SLs) N 
Wf(Pf t SLs) k (^ ~ Pf,SLs) N k ]■ 


k , 


(26) 


Thus, the optimal fusion rule, k opt , which minimizes the Bayesian risk can be obtained as 


h opt 


In 


'W„ 


1 -P‘ 


f,SLS 


N 


Wf ‘ l ~ P ' d ,SLS 


) 


ln( 


^ P d,SLS^^ P 'f,SLS ) 

T P d,SLS^ ( P 'f,SLS ) 


(27) 


where |~.~| denotes the ceiling function. 

Fig. 3 shows the TER as a function of the normalized threshold A n over NL and RL composite 
fading channels using different fusion rules. It is clear from the graph that the optimal fusion 
rule yields the best performance for both NL and RL channels. The minimum TER of the 
optimal fusion rule (k = 3) in the NL channel is 10 -3 9 , whereas in the RL channel it is 
lower by an order of magnitude which is 10 -2 9 . On the other hand, the AND rule yields the 
worst performance, yielding a minimum TER of 10 -1 for the RL scenario and 1(L 2 for the 
NL scenario. Interestingly, the OR rule did not improve from RL to NL and the minimum TER 
islCL 2 for both channels. Therefore, we conclude that using optimal fusion rule always increases 
the performance as compared to other fusion rules. 

In Fig. 4, we employ the same fading scenario with having different weights for miss-detection 
and false alarm, M = 2 and 70 = 5 dB. As observed, the optimal fusion rule holds its optimality 
with a minimum TER of 10 -4 5 for the NL channel and 10~ 3 ' 8 for the RL channel. The OR rule 
improved slightly with a minimum TER of 1CT 2 ' 8 . Finally, the performance of AND rule is still 
the worst amongst all the fusion rules in both channels. 
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Fig. 3: Total error rate vs normalized detection threshold over the NL channel using different 
fusion rules with 7 0 = 10 dB , m = 4, ( = 0.5 dB, M — 1, N — 10, W m = 0.3, Wf = 0.7 , 

and q = 10 -2 . 


B. Optimal Number of Antennas 

In this section, we investigate the effect of increasing the number of antennas with SLS 
diversity reception. It has been noted that SLS can increase the performance of the CSS remark¬ 
ably. Here, we found out that the performance is not monotonically increasing if the feedback 
channel is erroneous. Fig. 5 shows the effect of the feedback channel on the TER while increasing 
number of antennas M in various composite fading channels. In the perfect feedback channel 
case, the TER is monotonically decreasing as we increase the number of antennas. However, in 
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Fig. 4: Total error rate vs normalized detection threshold over the RL channel using different 
fusion rules with 70 = 5 dB, ( = 0.5 dB, M = 2, N = 10, W m = 0.7, Wj = 0.3 , and 

q = 10 " 2 . 


the imperfect feedback channel case, the total error rate decrease until a specific lower bound 
whereby increasing M afterwards is ineffective and inefficient. Since the function in the imperfect 
feedback channel case is convex, the optimal number of antennas M* can be evaluated when 


dQ m dQ f 
dM + dM 


(28) 
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Fig. 5: Total error rate vs number of antennas over perfect and imperfect feedback channels 
with 7o = 5 dB, C = 0.5 dB, k = 3, N = 10, A n = 8 , W m = 0.3, W f = 0.7 , and q = 10~ 2 . 


where 


DRAFT 


dQf 

dM 


N 

E 

n=k 


( N ^ 


n(P' ) n-1 f’ SLS /-I _ p> \ 
n Vf,SLS) Q M l 1 rf'SLS) 


N-n 


\ n / 

N ( N ^ 


E 

n=k 

dP' 


\ n 


(P'f,SLs) n (N ~ n)(l — P'f^SLs) 


TV—n—1 


X 


f,SLS 


dM 


dP\ 


f,SLS 


N f N \ 


dM 


•E 

n=k 


( P'flSLsY (1 — P'f,SLs) 


N-n 


\ U 


(n — (N — n). 


P\ 


f,SLS 


1 ~P\ 


f,SLS 
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Fig. 6: Optimal number of antennas vs. average SNR in various fading channels with, k = 3, 
N = 10, A n = 8, W m = 0.7, W f = 0.3 , and q = 10" 2 . 


and 


dP'f,SLS 

dM 


(l-P / ) M (-l + 2g)ln(l-P / ). 


(30) 
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Similarly 


dQm 

dM 


N ( N \ 


E 

n=k 


n(P: 


d,SLS 


r)P' 

,n— 1 d,SLS f -j 

1 dM 1 d ’ SLs) 


\ H / 
N ( N \ 


E 

n=k 


n y 


N - n)(l - P’t SLS ) 


N—n—1 


dPr 


X 


d,SLS 


dM 


dP’ 


d,SLS 


N { N \ 


dM 


E 

n=k 


iPU.SLs) n - 1 (l-P' d ,SLs) 


N-n 


\ / 


X 


(n — (N — ri). 


P' 


d,SLS 


1-P 


d,SLS 


and 


(31) 


—J\J— “ ( 1 _ Pd,Mc) M q ln(l - Pd, mg) 

-(l-P d , MG ) M (l-q) ln(l — Pf) 

= (1 — P d ,M G ) M (—1 + 2g) ln(l — Pd, mg)- (32) 

The optimal number of antennas can be obtained by substituting (29) and (31) in (28) and solving 
it for M numerically. 

Fig. 6 shows the optimal number of antennas as function of the avarage SNR, over various 
composite and non-composite fading channels. As expected, the optimal number of antennas is 
inversely proportional to 70 . The RL channel requires in general more antennas than the NL 
and Weibull channels to reach M*, since the NL and Weibull channels depicts mild multipath 
fading with nn = 4. 


DRAFT 


March 4, 2015 
















ALHAMMADI et al : UNIFIED ANALYSIS OF COOPERATIVE SPECTRUM SENSING OVER COMPOSITE AND GENERALIZED FADING CHANNELS 19 


VI. Conclusions 

We have investigated the performance of cooperative spectrum sensing over generalized fading 
channels. A unified MG based approach is proposed for the performance analysis of cooperative 
spectrum sensing in cognitive radio. It has been found that in a composite fading channel, the 
performance can be improved remarkably by applying an optimal fusion rule outperforming 
the AND rule. Moreover, increasing the number of antennas in a CSS network with erroneous 
feedback channel does not decrease the TER monotonically. Based on this observation, the 
optimal rule for number of antennas is determined numerically. Finally, it has been found that 
the fading severity can highly affect the optimal rule for the number of antennas. 
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